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We explain an accelerated expansion of the present universe, suggested from observations of 
supernovae of type la at high redshift, by introducing an anti-frictional force that is self-consistently 
exerted on the particles of the cosmic substratum. Cosmic anti- friction, which is intimately related to 
"particle production" , is shown to give rise to an effective negative pressure of the cosmic medium. 
While other explanations for an accelerated expansion (cosmological constant, quintessence) in- 
troduce a component of dark energy besides "standard" cold dark matter (CDM) we resort to a 
r/-\ ' phenomenological one-component model of CDM with internal self-interactions. We demonstrate 

l/~) , how the dynamics of the ACDM model may be recovered as a special case of cosmic anti-friction. We 

(f) ■ discuss the connection with two-component models and obtain an attractor behavior for the ratio 

of the energy densities of both components which provides a possible phenomenological solution to 
^^ , the coincidence problem. 
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I. INTRODUCTION 



There is evidence from SN la data for our present universe to be in a state of accelerated expansion ]lHq|. This 
interpretation, which is indirectly also backed up by recent data from the balloon experiments Boomerang |6J and 
Maxima M according to which we live in a flat universe, requires a cosmic medium with sufficiently high negative 
pressure to violate the strong energy condition (SEC), p + 3P > 0. Cosmic matter with negative pressure is now 
known as "dark energy" . The problem why the density of the dark energy is of the order of the matter density just at 
the present epoch is the "coincidence problem" [^0| . An obvious dark energy candidate is a cosmological constant of 
the order of the current critical density. Another option is a scalar field called "quintessence" , either with a suitable 
potential [p|-p7[ or with a nonstandard kinetic term p8[ . 

Negative pressures also occur in a different context. As was first pointed out by Zel'dovich pit] and Hu [Bof, quan- 
tum processes in the early universe such as cosmological particle production may phenomenologically be equivalent 
to effective negative bulk pressures. Numerous investigations have subsequently explored this analogy |2l|-|29(|. For 
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gaseous matter with specific internal self-interactions negative cosmic bulk pressures can be derived within the frame- 
work of relativistic gas dynamics |3Q-|33|]. In particular, it turned out that accelerated expansion (inflation) of the 
universe can be driven by a self- interacting gas |p0|-p3|. 

The observational data of the early 90th left ample room for a presently non-negligible dissipative bulk stress on 
cosmological scales, as has been pointed out by two of us in Ref. p4j] , Hypothetically, this bulk stress was ascribed to 
internal interactions inside the dark matter, the latter being supposed the dynamically dominating component of the 
universe. According to the state of knowledge at that time our studies were restricted to matter for which the SEC 
holds. Taking into account the recent development sketched above, it seems natural to investigate a present phase of 
accelerated expansion from a gas dynamical point of view. 

The present work is based on the assumption that the observational evidence for an effective cosmological constant 
is an indication for the existence of additional interactions within the cosmic medium, which macroscopically manifest 
themselves as negative pressures. Our strategy is the following: Instead of introducing from the outset a new kind 
of matter ("quintessence" or "Q matter") with a negative pressure in addition to cold dark matter (CDM), we start 
with a one-component description of the cosmic medium. Although the corresponding substratum is assumed to 
consist of non-relativistic particles, it does not represent simple dust since we include interactions within the matter. 
These interactions turn out to be equivalent to effective one-particle forces, which are self-consistently exerted by the 
cosmic medium on each of its individual particles. We show that the cosmological principle restricts these forces to act 
like "friction" or "anti-friction" . Generalized equilibrium requirements for the cosmic medium relate those forces to 
temperature and chemical potential. For a non-relativistic substratum only cosmic anti-friction generates a negative 
fluid bulk pressure. We demonstrate that a suitable amount of cosmic anti-friction leads to a SEC violation equivalent 
to an accelerated expansion of the universe. On this basis we argue that the hypothesis of cosmic anti-friction offers 
an alternative interpretation of the magnitude-redshift relation for type la supernovae. 

We present three different phenomenological models of cosmic anti-friction which give rise to Hubble diagrams 
consistent with the SN la data. We estimate the redshifts at which the accelerated expansion started. Although SN 
la data alone cannot discriminate between those models, observations of the cosmic microwave radiation anisotropics 
rule out one of them. 

In a subsequent step we demonstrate how the anti-friction dynamics may be decomposed into a two-component 
picture with one component being Q matter, the other one CDM. This procedure reveals that a negative fluid pressure 
due to anti- friction may be dynamically equivalent to the accelerating effect of dominating Q matter. For a specific 
ansatz for the cosmic force and for a specific decomposition one recovers the ACDM model. Moreover, it turns out that 
the same amount of accelerated expansion of the presently observed universe is compatible with different splittings of 
the total energy density into the energy densities of Q matter and CDM. In particular, there are decompositions in 
which Q matter decays into CDM. For a specific decay rate there exists a stable attractor solution with a fixed ratio 
of the energy densities of CDM and Q matter, which indicates a possible solution of the coincidence problem. 

The paper is organized as follows. Section II recalls some basic relations and features on the possible role of a 
cosmic bulk pressure. In Sec. Ill we we provide the basic relations of a kinetic theory for self-interacting gases. We 
show that a certain class of particle number non-preserving interactions may be mapped onto effective one-particle 
forces. The cosmological principle implies, that these forces necessarily describe a friction or an anti-friction within 
the cosmic medium. We determine the strength of the force and realize that only anti-friction is compatible with a 
negative pressure. The back reaction of the anti-frictional self-interactions on the cosmological dynamics is considered 
in Sec. IV for three different anti-friction models. The luminosity and angular distances as functions of redshift are 
studied for all three models. Corresponding two-fluid models are established in Sec. V which provide the basis for a 
discussion of the coincidence problem. Sec. VI sums up our conclusions on the possible role of cosmic anti-friction. 
Units have been chosen so that c = ks = Ti = 1. 

II. COSMIC BULK PRESSURE 

We suppose the Universe to be describable by the stress-energy tensor of an imperfect fluid 

T$ = {p + P)u a u b + Pg ab , (1) 

where p is the energy density, measured by an observer comoving with the fluid 4-velocity u a normalized according 
to u a u a = —I. The effective pressure P splits into two parts, 

P=p+U, (2) 

p being the equilibrium pressure with p > for gaseous matter and II a non-equilibrium part. For a perfect fluid we 
have II = 0, i.e. P > 0. For a conventional viscous fluid II < is valid during expansion, e.g., if kinetic energy of the 



fluid is transferee! to internal degrees of freedom. In the first-order Eckart theory one has II = —3H( (see, e.g. p5[), 
where H = u a a /3 = a/a is the Hubble expansion rate and £ > is the coefficient of bulk viscosity. Within the more 
satisfactory second-order theories II becomes a dynamical degree of freedom (see, e.g., [p6|~p9[). Both the first- and 
the second-order theories are valid under the condition | IX | < p, such that the effective pressure P of a viscous fluid 
or gas is positive. 

Apart from viscosity, particle number non-conserving interactions inside the matter may, as discussed above, lead 
to an effective bulk pressure. This includes particle production out of the gravitational field. The fact that II < if 
there is particle production and that P < is possible under such conditions may be demonstrated as follows. Let 
the cosmic matter be characterized by the particle flow vector 

N l = nu l , (3) 

where n is the particle number density. In case the fluid particle number is not preserved, the number density changes 
according to the balance p3| , p4| 

N? a = h + 3Hn = nT , (4) 

where T = N/N is the change rate of the number N = na 3 of particles in a comoving volume a 3 . For V > we have 
particle creation, for T < particles are annihilated. Conservation of the effective stress-energy tensor (fy) implies 

p + 3H{p + P) = 0. (5) 

With the help of the Gibbs equation 

Tds = d^-+pd-, (6) 

n n 

where s is the entropy per particle and T the temperature, and with the balances (H) and (ra) we obtain 

nTs = -3HU -{p + p)T . (7) 

If the particle number N is conserved, i.e. for r = 0, the second law of thermodynamics implies II < in an expanding 
universe. 

If the particle number is not conserved one may define "isentropic" (or "adiabatic") particle production by s = 0, 
which here means "constant entropy per particle" . Under this condition the equilibrium entropy per particle does 
not change as it does in dissipative processes. Instead, one can associate a viscous pressure to the particle production 
rate ||Q: 

s = =► n = -(p + p)^. (8) 

The cosmic substratum is not a conventional dissipative fluid but a perfect fluid with varying particle number. 
Obviously, T > guarantees II < 0. Substantial particle production is a phenomenon which is reasonably to be 
expected in the early universe. It is less clear whether such processes are operative at the present epoch as well. 
However, given that the nature of CDM is unknown, there seems to be some room for speculations in this direction. 
As easily seen, P < in the case of dust (p — 0) and P ~ — p is possible if T/(3H) = 0(1). Compared to typical 
rates of particle physics this requires an extremely small particle production rate only. Below we shall comment on 
the origin of negative bulk pressure in more detail. 

Referring to matter creation as a relevant cosmological mechanism may remind of corresponding processes within 
the steady state model [Eoi. However, different from the latter, our considerations are fully within Einstein's theory. 
Moreover, we shall trace the production process to internal interactions within the system. 

For a universe of the Friedmann-Lemaitre- Robertson- Walker (FLRW) type with scale factor a we have 

2-+(-) +A = _ 8 ^ G p, (j. = 1,0,-1) (9) 

a \a J a z 
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p 8nG 

= n-i, n^^- = — Pl (io) 



a 2 H 2 ' p r 3H 2 



implying 



n = - 



4ttG 



1 



Q - q )H 2 + 4irGp 



(11) 



where q = —a/(aH 2 ) is the deceleration parameter. In the standard big-bang scenario the non-equilibrium pressure 
II is ignored and since p cannot become negative it folllows £l(t) < 2q(t). But it is obvious that as long as p(t) does 
not vanish and both pressures are of similar magnitude, fl(t) can be either larger or smaller than 2q(t) . Usually the 
current value of the hydrostatic pressure is approximated by the state equation of dust Pq — 0. The recent SN la 
data seem to confirm the existence of an effectively negative pressure of the cosmic medium. Assuming pq ~ in Eq. 
(O) , the current value of the bulk stress can then be expressed as 
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where the current energy density po is of the order of the critical energy density p c o = 3Hq / (8irG) . For go > we get 
3|IIo| < po, which is in agreement with the strong energy condition. The latter is violated, however, for go < 0- As 
already mentioned, for "conventional" viscous matter without particle production the non-equilibrium part II of the 
pressure is smaller in magnitude than the equilibrium contribution p, so that P > 0. Here we argue that relation m) 
offers the option, to understand the existence of a negative bulk pressure as a manifestation of cosmological "particle 
production" . In the following sections we show how such a kind of negative pressure may emerge as a consequence of 
specific internal interactions within the cosmic substratum. 



III. KINETIC THEORY FOR SELF-INTERACTING GASES 



A. Basic relations 



The one-particle distribution function / = / (x,p) of a relativistic gas is supposed to obey the Boltzmann equation 

(13) 



Jc / df 



L [/] = p*f,i ~niP k P l g^ = C[f]+S (x,p) , 

where C[f] is Boltzmann's collision integral. The term S(x,p) on the right-hand side takes into account additional 
interactions which can not be reduced to elastic, binary interactions. In particular, it describes production or decay 
processes of particles. 

The particle number flow 4-vector N l and the energy momentum tensor T lk are defined in a standard way (see, 
e.g., gTJ) as 



N' 



J dPtff (x,p) , T lk = J &Pp l p k J (x,p) 



(14) 



The integrals in the definitions (|4|) and in the following are integrals over the entire mass shell, characterized by 



VVi 



and p° > 0. The entropy flow vector S a is given by 41 



S a = -JdPp a [f\nf-f] 



(15) 



where we have restricted ourselves to the case of classical Maxwcll-Boltzmann particles. 
Using the general relationship |f3[ 



p ai ....p a "p b fdP 



Jp a K..p a "L[f]dP 



and eq.(|13|) we find 



N^ a = JdP(C [/] + S) , Tf k = J dP P a (C [/] + S) , 



(16) 



(17) 



and 



S°< a = a c + a s , (18) 

where 

a c = -JdPC[f]lnf, (19) 

and 

a S = - fdPSlnf. (20) 



Under the condition that with respect to the elastic part of the interactions the gas is at equilibrium, the expres- 
sion In/ is a linear combination of the collision invariants 1 and p a and the collision integral C [/] vanishes. The 
corresponding equilibrium distribution function becomes (see, e.g., |4lf| ) 

f°(x,p)=exp{a + f3 aP a } , (21) 

where a = a (x) and B a (x) is timelike. This implies ac = 0, i.e., there is entropy production only due to as- 

With / replaced by /° in the definitions (H) and (113), N a , T ab and S a may be split with respect to the unique 
4-velocity u a according to 

N a = nu a , T ab = pu a u b + ph ab , S a =nsu a . (22) 

Note that we have identified here the general fluid quantities n, p and p of the previous sections with those emerging 
from the Maxwell-Boltzmann gas dynamics. The exact integral expressions for n, p, and p may be found, e.g., in Ref. 
El . The entropy per particle s is 

,= e±*-±. (23) 

nT T y ' 

Here we have used the identifications /3j = 0Ui, 8 — T _1 , and a — plT with p being the chemical potential. 
Use of the equilibrium distribution function (Elf) in the balances (117]) yields 



n + 3Hn = nT = / dPS" , (24) 

and 

u a [p + 3H(p + p)} + {p + p)ii a + p ib h ah = -t a = f dPp a S° , (25) 

where S° denotes the source term S for / = f°. We consider the special case that 5° depends linearly on f°. For 
reasons that will become clear shortly, we suppose that the factor of proportionality can be written in terms of a 
suitable projection of a quantity F 1 which will turn out to play the role of an effective one-particle force: 

S° = -mdiFtf . (26) 

The constant factor — m has been chosen for later convenience. The expression (Eq) may be regarded as a special case 
of the more general structure 

S = -mF 1 ^- . (27) 

dp 1 K ' 

It is straightforward to realize that a "collision" term of this form may be taken to the left-hand side of Boltzmann's 
equation ([13|), resulting in 

P7,i-nyyJ£+mi^ = C7[/] . (28) 

The left-hand side of this equation can be regarded as 

df(x,p) = df dx l df dp 1 
dX dx l d\ dp 1 d\ 



with 

Equations (£9|) are the equations of motion for gas particles which move under the influence of a force field F l — 
F l (x,p). The quantity A is a parameter along the particle worldline which for massive particles may be related to the 
proper time r by A = t/to. Consequently, a specific "collisional" interaction, described by a "source" term S, may 
be mapped onto an effective one-particle force F l . This demonstrates that there exists a certain freedom to interpret 
collisional events in terms of forces. (This freedom can also been used in the reverse direction, i.e., to interprete (parts 
of forces) as collisions |4q| ) . We emphasize that our approach is different from the "canonical" theory of particles in a 
force field for which the force term mF l df /dp 1 in Eq. fll9| ) is replaced by md (i 714 /) /dp 1 |45| . While both approaches 
are consistent with the equations of motion (J29J), they coincide only for OF 1 /dp 1 — 0, which holds, e.g., for the Lorentz 
force. In the cases of interest here we will have OF* /dp 1 ^ 0. 

B. Cosmic forces 

Since we assume the universe to be homogeneous and isotropic at large scales we ask for forces which are consistent 
with the cosmological principle. In such a case the metric of space-time is of the Robertson- Walker form and the energy- 
momentum tensor is given by Eq. (HI). Since the particle four-momenta are normalized according to p % pi = — m 2 , the 
force F l has to satisfy the relation piF % = 0. The momentum of a comoving particle is pi-. — mu 1 . From its definition 

comoving particles are force- free and thus F l (mu) = 0. This property follows also from the relation 

^ =«?«-. (30) 

dr ' to 

which for pi-. = mu 1 via u l . n u n = u 1 oc F % requires a vanishing force since the cosmological principle implies u l = 0. 

On a spatial slice T, t , normal to u l , the force field F l has to be independent of the spatial position, otherwise 
homogeneity would be violated, thus F l — F l (p;t). Since p l and u 1 are generally independent 4- vectors, we may 
decompose according to 

F* = Ap l + Bmu 1 , (31) 

where A and B are arbitrary functions of p 1 with dimension 1/time. From F l pi = we find 

F l = — {-Ep l + m 2 u l ) , (32) 

to 

with E = —p l Ui being the particle energy as measured by a comoving observer. For a comoving particle one has 
E = m and we consistently recover that F l (mu) = for all B. A particle which exactly moves with the mean 
macroscopic four-velocity is force free. The temporal and spatial projections of the force field give 

UiF* = —(E 2 - m 2 ) , e i F i = -—E^E 2 -m 2 , (33) 

m m 



where 



Ve 2 



(p l - Eij) (34) 



is the spatial direction of the particle momentum (e l u.i — 0, e'e,; = 1). Due to spatial isotropy B may not depend 
on the spatial direction e 4 , thus B = B(E; t). The expression ( |32j ) is the most general force field consistent with the 
cosmological principle. 

According to the projections (B3h, the force is acting parallel or anti-parallel to the motion of the particle under 
consideration, depending on the sign of B. For non-relativistic particles it should be a good approximation to assume 
that B is independent of E. With E = m + e where e = rav 2 /2 <C to, we find at leading order in the velocity 
eiF 1 ~ —B(m)mv. This is nothing but Stokes' law of friction. For B > the force field may be interpreted as cosmic 
friction, for B < as cosmic anti-friction. We conclude that cosmic (anti-)friction is the most general force field which 
is compatible with the cosmological principle. 



With the expressions ( p6| ) for <S° and the equivalent force ([32|) we may calculate the "source" terms in the balances 
( J24| ) and (|2^). For a B independent of E, the results are 

T = -2.B , (35) 

and 

t a = 3Bu a {p + p) . (36) 

In general, neither the particle number nor the energy momentum are conserved. Consequently, T ab in Eq. ( [2 4 ) is 
not the quantity which will appear on the right-hand side of Einsteins' field equations. Eq. (pq) clarifies that the 
force strength determines the particle production (decay) rate T = N/N. The entropy production density (|1S|) is 
determined by 

S?a = °"5 = -aN? a - (3 a T ab b = nsT = -3nsB . (37) 

A production of particles is characterized by T > and corresponds to B < 0, i.e., an anti-frictional force, while a 
decay of particles is equivalent to an effective friction. 

C. Generalized equilibrium solutions 

To obtain the conditions under which the equilibrium distribution (|2lj) is preserved even under the action of an 
(anti-)frictional force we insert the expression (plf) into the Boltzmann equation (Eq) which yields 



P a a, a + P(a;b)P a p b = -raftF* . (38) 

Since /3j = Ui/T, it is only the projection uiF 1 of the force which is relevant here. If this projection vanishes, relation 
(pq) reduces to the "global" equilibrium condition of standard relativistic kinetic theory, i.e., to a = const and either 
to the Killing vector condition f3t a -b) — for m > 0, or to the conformal Killing vector condition /3( a; M = 4> i x ) 9ab 
for m = 0. The Friedmann models do not supply a timelike Killing vector, thus there is no equilibrium solution 
for particles with finite mass. However, in the non- relativistic limit T ^C m a quasi-equilibrium, characterized by 
a = m/T + const, Pi = Ui/T and T ex a~ 2 exists. In the case of cosmic (anti-)friction we find a similar quasi- 
equilibrium solution for Friedmann models which reads 
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Bj . (39) 

With 3B = — r = N/N according to relation (p5|), the temperature behavior is 

T ex a- 2 N 2 / 3 . (40) 

For vanishing (anti-)friction the particle number is constant and the familiar T ex a~ 2 dependence for nonrelativistic 
matter is recovered. 

The explicit knowledge of the force F l in terms of Y now allows us to study the motion of the matter particles 
explicitly. Contracting the equation of motion (E9) with the macroscopic four-velocity results in 

T)(u lP l ) dE _, 1 i k 

-^ =—r-= UiF 1 + —Ui- k p l p k . 41 

dr dr m 

With Eq. ( p3| ) and under the condition of spatial homogeneity we have 

dE B + H 



dr 



(E 2 -m 2 ) . (42) 



Since dr = dt (m/E) and dE/dt = E where E = m + e with e ^C m, we finally obtain for the evolution of the 
non-relativistic (kinetic) energy the expression, 

e = -2 (h - D e . (43) 



Consequently, 



e ex a~ 2 N 2/3 ex T , 



(44) 



the equipartition theorem. Since the exponent a + f3 a p a = a — E/T of the function (21) in the nonrelativistic limit, 
with a = m/T + const and E = m + e, reduces to a — E/T — > const + e/T, relation (44) demonstrates explicitly the 



invariance of the equilibrium distribution (plj). The nonrelativistic velocity scales as v 
With Eqs. ( p5| ) and ( |3q ) the fluid energy balance becomes 

p + 3H(p + p) = -3B(p + p)=T(p + p) 



tx a 



LjV 1 / 3 . 



Introducing the quantity 



U=-(p + p) = -—(p + p) 



(45) 



(46) 



Eq. ( flq ) may be written in the form of the energy balance (ra), where II plays the role of an effective viscous pressure 
according to the definition (|2j), thus we have P = p + H. The energy balance for T ab with generally non- vanishing 
source in ( pq ) is identical to the conservation law u a T^. k = with the energy momentum tensor p). It is this energy- 
momentum tensor which appears on the right-hand side of Einstein's field equations (see the discussion following Eq. 
(pq)). With the re- interpretation (Bfl) of the sources on the right-hand side of the balance (E3) in terms of an effective 
pressure of the medium the latter becomes a "closed" system. In the following we shall restrict ourselves to negative 
pressures II, corresponding to the production of particles and to an anti-frictional force. 

We emphasize again that despite of the non- vanishing entropy production S a a the microscopic particles are always 
governed by an equilibrium distribution function. In this context S a a > describes just an enlargement of the phase 
space of the system but not a dissipative process. Although inter-particle collisions are necessary to establish an 
initial equilibrium characterized by Eq. ( pl| ) , this equilibrium may then be maintained even under the influence of the 
anti-frictional force and in the absence of further collisional interactions which might have been frozen out. The force 
which gives rise to a negative fluid pressure is compatible with an equilibrium distribution of the particles during the 
expansion. This feature is an essential advantage of the presented approach since it allows us to apply standard gas 
dynamical concepts to characterize "exotic" matter forms. 



IV. COSMOLOGICAL DYNAMICS 



Let us assume the cosmic substratum after matter-radiation decoupling to be non-relativistic matter with internal 
anti-friction, characterized by the energy-momentum tensor (fil) with p <C p and 



3H H H H 

Bearing in mind that T = N/N the energy balance (M5T) may be integrated to yield 

N /a \3 



P = Po 



N \a) 



(47) 



(48) 



The index again denotes the present epoch. According to the Friedmann equation for the spatially flat case, 
8irGp = 3H 2 , to which we restrict ourselves from now on, the corresponding Hubble rate is given by 



H = H 



N /ao^ 3 
N Q \a 



1/2 



(49) 



For the ratio H/p we obtain 



n _ b _ b 



No ( a 



N \a 



1/2 



(50) 



The "particle number" changes as 



N = N Q exp 



-3 



dt\B\ 



(51) 



For vanishing anti- friction, corresponding to a conserved particle number, we have N = N . 
The ratio \B\/H enters the Hubble law for small redshift z (dj, is the luminosity distance), 

H d L = z + - (1 - q ) z 2 + .... , 

via the deceleration parameter q, which for k = and nonrelativistic matter becomes 

1_3 M _ 
1 2 2 H 

Generally, the luminosity distance d^ in a spatially flat universe may be written as 

d L = (l + z) —— 
Jo H{z) 

As usual [|46| this is related to the angular distance d\ of an object as 

d A (z) = (l + z)- 2 d L {z) . 

Relation (f47|), the field equations @ and (JHJ) for k = and p <C p may be combined to yield 

n s 2fl 

7 ~ h ~ ~ ~ HJp ' 

It is convenient to express the Hubble rate as a function of redshift z = (ao/a) — 1. With 

H = -H'H(1 + z) , 
where H 1 = dH/dz, the resulting equation is 

H' 3 

B + H ~ 2(1 + z) ' 



(52) 



(53) 



(54) 



(55) 



(56) 



(57) 



To establish specific models, assumptions about B are necessary. In the following we solve the cosmological dynamics 
for three different choices of B and compare the results with the Hubble diagrams from recent SN la measurements. 
Although SN la data cannot discriminate between our different models, since they provide equally good fits for redshifts 
that are accessible with SN of type la, the additional information from the angular distance to the last scattering 
surface and CMB measurements will rule out one of the choices even without any detailed statistical analysis. 

A. The case B ex -H 

The simplest possible choice is apparently B = —vH with v = const. According to Eq. ( |47| ) this ansatz is equivalent 
to assuming a constant ratio H/p, which should be considered realistic at most piecewise. With this assumption we 
get a Hubble rate 



H(z) = H (l + z)^ 1 -^ 



(58) 



Note that the same power-law behavior follows if a perfect fluid with the equation of state —v = w= p/p is assumed. 
We recover the matter dominated universe for v = and the vacuum dominated universe for v = 1. Equation (J5J) is 
easily integrated to provide an explicit expression for the luminosity distance (y ^ 1/3) 



H Q d L 



3u- 1 



(1 + 2) 



1(1+3*) 



1-2 



(59) 



This result coincides with the one by Lima and Alcaniz |47| . The corresponding luminosity distance - redshift relation 
is shown in Fig. |l] while Fig. g shows the difference to the Standard Cold Dark Matter (SCDM) model. We find that 



v = 0.5 gives a good fit to the observations of SN la. At first sight one might have expected the best fit for v = 0.7 
since it is this value which reproduces an equation of state P/p = -Po/po = —0.7. However, as already indicated, 
a constant ratio P/p is realistic only piecewise. The fact that the better fit is v = 0.5 suggests that v must have 
been smaller for z > than for z = 0. This illustrates the circumstance that the SN la observations for higher z do 
not directly reflect the cosmological equation of state at the present time but, as to be seen from the integral in the 
expression ( p4j ) for dt, depend on the entire dynamics from redshift z to redshift 0. 
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FIG. 1. Hubble diagram for (SIm^a) = (1,0) [thin, dashed line], (0.3,0.7) [thin line] compared to case A with the values 
v — 0.7, 0.5, 0.3 [thick lines from top to bottom]. The data points are taken from Perlmutter et al. 1999, the diamonds are the 



Calan/Tolodo SNla data, the triangles are those of the SCP. 
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FIG. 2. Differences of the magnitudes with respect to a (Qm, J7a) = (1, 0) universe versus redshift. The plotted models and 
data are the same as in figure 1. Note that other authors often depict the magnitude difference with respect to an empty 
universe, which has negative spatial curvature and is thus incompatible with inflation. 

The "particle number" changes according to 



N 



exp 



|B| dz 
"ff (1 + z) 



(1 



-3v 



(60) 



The energy density [cf. Eq. (p|)] decays as p ex (1 + z) 3 '- 1 ^^, the temperature [cf. Eq. (|(i|)] as T oc (1 + z) 2{ - 1 v \ 
For the 'best fit' v = 0.5 there is a decrease in N/N by about 3(5) for a redshift of 1(2). Over a larger interval N 
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might decrease by many orders of magnitude. This again reflects the unrealistic nature of F '/ ' p = const over a large 
range of z values. 

The observations of SN type la are restricted to redshifts of order 1. In order to test a specific model of anti- friction 
at redshift z»l, additional information is required. A very promising possibility are CMB anisotropy data. Acoustic 
oscillations in the early universe give rise to peaks and dips in the band power spectrum of these anisotropies. The 
typical length scale of the largest acoustic oscillations is given by the sound horizon R s = c s /H at the time of photon 
decoupling, where c s is the sound velocity. This physical length corresponds to the presently observed angular scale 
pq| of the first acoustic peak 

RsjZdcc) ,„-.-. 

Olst peak ~ ~r~, T • (01,) 

«A(2decJ 

In fact, the corresponding observational data rule out any constant value for P/ p. Namely, insertion of Eqs. (pq), 
( J5q ) and ( |59| ) leads to (zdcc ^> 1) 

<*ist peak ~ c s — - — « 0.14rad , (62) 



which for v = 0.5 and c s = y 1/3 exceeds the detected <5i st pC ak — 0.9° = 0.016 rad [|6||7]] by about one order of 
magnitude. Thus, a constant effective equation of state over a large range in redshift is incompatible with observations. 
To illustrate this feature we plot in Fig. ffl the angular scale under which the Hubble radius H^ 1 (z) is seen for the 
'best fit' v = 0.5. At z — Zdcc ~ 1100 the Hubble radius differs from the sound horizon only by a factor c s « yi/3. 
The present model corresponds to the upper curve. In the following we consider two cases for which P/p is time 
varying. Since we are interested in a modification of the cosmological dynamics at late times we focus on two models 
where this ratio is increasing with time. 

B. The case B — — \Bq\ = const 

From Eq. fl47|) follows that the simplest choice leading to an increasing ratio |Hj/p is a constant value of B. From 
Eq. d57J ) we obtain with v = \B \/H , 

H(z)=H [(l-p )(l + zf' 2 + uo] ■ (63) 

This solution describes a transition from matter domination to vacuum domination of the universe. In the past, z ^> 1, 
we have H oc (1 + z) 2 ' 3 , a matter dominated universe, while for the future, z — ► — 1, H and thus p are constant. The 
expansion of the universe starts to accelerate when H(z) < 3\B\ [cf. Eq.(p3)], which happens for a redshift z acC j given 
by 



l + z a 



2^o 

1 - ISO 



2/ .3 

(64) 



In Figs, pfl and we compare the Hubble diagram for this model to recent observations of SN la. Here we find a 'best 
fit' value vq = 0.7. For this model the universe starts acceleration at z acc sa 1.8. The increase in the "particle number" 
may be expressed as 



i2 



No (l + z ) 

which gives an increase by a factor of about 5 since z acc . 



N [(l-^o)(l + z) 3 / 2 + ^o] 

— " 3 ' ( 65 ) 
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FIG. 3. Same as in figure 1 but for case B with \Bq\/Hq = 0.9, 0.7, 0.5 [thick lines from top to bottom]. 
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FIG. 4. Differences of the magnitudes versus redshift for the models in figure 3. 

Also for this model the angular scale of the Hubble radius is shown in Fig. for i>q = 0.7 (second curve). Obviously, 
here the Hubble scale is of the order of the scale of the first acoustic peak, i.e., this model is consistent with the CMB 
observations. However, a more quantitative statement could only be made on the basis of a Boltzmann code, which 
is beyond the scope of our present work. 

C. The case B ex -H' 1 
Another choice leading to an increasing ratio |n|//o is \B\ ex H , equivalent to an ansatz 



IS! 



1 Hi 



H n + 1 H 2 ' 

The specific choice of the constant factor 1/ (/i + 1) was made for later convenience. Integration of Eq. ([ 
ansatz (|66|) yields 



(66) 
with the 



H=j^= w+zf+ ir 



(67) 
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For z>lwe have 



H oc (1 + z) 



3/2 



(68) 



which is characteristic of a matter-dominated universe. For the opposite case z — * — 1 the Hubble rate approaches the 
constant value 



ff 



H n 



V7h~T 



(69) 



The Hubble rate (|67]) implies again a transition from a matter-dominated universe at z > 1 to a de Sitter universe 
as z — > — 1. Equations ( p7|) and (54) determine the luminosity distance, which is plotted in Figs. || and g. This case 
in fact includes all ACDM models, as can be easily seen by replacing l/(/i + 1) — ► SI a and /x/(/x + 1) — > ^cdm- As 
expected, the 'best fit' model has 1/(1 + /i) = 0.7. 

To obtain the redshift at which acceleration starts we write 



H 



/i(l + z) 3 + l 



(70) 



The condition \B\/H > 1/3 for accelerated expansion [cf. Eq. (|5^)] gives 1 + z acc = (2/fi) 1 ^ 3 . The growth of N 
follows from 



N = N, 



^(l + z) 3 + l 

(/i + l)(l + z)3 



(71) 



For the favored value 1/[1 + /i) = 0.7 accelerated expansion starts at z acc s=s 0.67, which is in the expected range 
0.5 < z acc < 1 [cf. Ref. p8[], and A^ has grown since then by a factor of about 3. 
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FIG. 5. Same as in figure 1 but for case C with 1/(1 + /j.) — 0.9, 0.7, 0.5 [thick lines from top to bottom]. Note that this case 
is identical to ACDM, with £Ia = 1/(1 + fj,). 
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FIG. 6. Differences of the magnitudes versus redshift for the models in figure 5. 

Again we show an estimate of the angular scale of the Hubble radius (for 1/ (1 + /i) = 0.7) in comparison to the scale 
of the first acoustic peak in Fig. p] (lower curve) . This model is consistent with observations as well. The difference 
to the former model is large enough that it should give rise to a significant difference in the CMB predictions. We do 
not further elaborate this point here. 
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FIG. 7. The angular scale H~ 1 (z)/(1a(z) in degrees under which the Hubble radius H~ x is seen for the models of subsections 
IV. A., IV. B., and IV. C (from top to bottom). For comparison we indicate the angular scale of the observed first acoustic 
peak (~ 0.9°) by the horizontal line. 

V. TWO COMPONENT MODELS AND THE COINCIDENCE PROBLEM 

A. Interacting fluids 

In this section we investigate in which sense cosmological anti-friction may be regarded as effective description for an 
underlying two-component model. To this purpose we decompose the total energy density (^8() into a "conventional" 
matter part (subscript M) and a "Q matter" part (subscript Q): 



P = PM+PQ , 



(72) 



with 



14 



Pm = rp [ — 
V a 



and 



PQ 



N 



Po 



m 



(73) 



where r = r (a) > with r < N/N and r = r (a ) < 1. The factor r fixes the ratio of both components at a^: 

Pm (qq) _ r a 
PQ (oo) 1 - ro ' 



(74) 



For r — const the universe evolves as though the cosmic substratum would consist of two non- interacting components: 
Non-relativistic matter with pm — and Q-matter with an equation of state 



B 



\B\ N/N 
H (N/N ) - r 



Pq , (r — const) 



(75) 



For an expanding universe the pressure II is always negative. For any form of the potential of a quintessence model 
the corresponding coefficient of anti-friction is easily obtained from Eq. (F7q) . 
If r is allowed to vary we have 



Pm + SHpia = -pm ■ 



(76) 



For r > the term on the right-hand side may be regarded as a matter "source" . This corresponds to a "sink" in the 
energy balance of the Q component, 



PQ + "ZHpq =T(p Q + p M ) Pm ■ 



By introducing the effective pressures 



and 



Pm = n M 



r/r 



Pm 



p q = n - n M 

the balances ( |76| ) and (J77J) become 



T N r 



PQ 



(N/N ) - r ' 
p M + 3H {p M + Pm) = 



and 



p Q + 3H(p Q + P Q ) = Q, 



(77) 



(78) 



(79) 



(80) 



(81) 



respectively. The circumstance that for r > both components have a negative effective pressure (given that |IIj\/| < 
|II|) is essential for a discussion of the "coincidence problem" in the present context. 

Supernovae la observations suggest that \P\/p = \B\/H is of the order unity in the present epoch. Within our 
approach this raises the question: Why is the coefficient of antifriction —B of the order of the Hubble rate just at 
the present epoch? As demonstrated in Sec. IV. A. for the case B oc —H, the CMB data rule out that this was 
the case through the entire evolution of the universe since decoupling. Thus it seems to be a coincidence that we 
live in this special era. The decomposition (JT2|) allows us to relate our consideration to the usual discussion of the 
coincidence problem. It is obvious that all models according to which the universe is made of a non-interacting 
mixture of quintessence and CDM correspond to r = const. In this special case our approach does not provide any 
new insight into the coincidence problem. For r > 0, however, it offers a solution which is similar to the one proposed 
by Chimento et al. p9|. The quantity of interest is the ratio pm / PQ which is governed by the equation 



Pm 
PQ 



= 3H 



Pm 
PQ 



PQ 



Pm 
Pm 



(82) 



For Pm = 0, or equivalently r — const, one has Pq = II < 0, i.e., the ratio pm/pq continuously decreases and for 
large cosmological times one has pm <C pq. In other words, the matter component becomes dynamically negligible. 
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However, if an exchange between both components is admitted, which amounts to a nonvanishing quantity Pm — IIm, 
there exists a second stationary solution of Eq. (B3), namely 



Pq = Pm 

PQ PM 

Combining relations ( f78J ) and (J7J3) we obtain 

Pq Pm _ 1 N/N 



p Q p M SH (N/No) - r 
Obviously, the condition (B3h is equivalent to 



r 



r _ N 

r~N 



r = r 



N 
No' 



Via Eq. ( |73|) the stationarity condition (|89) for the ratio pm/pq provides us with 



Pm 

pq 



■'■<) 



l-r ' 



(83) 



(84) 



(85) 



(86) 



which, according to the split (|73|), is just the present ratio of both components. Concerning its possible role for the 
coincidence problem it is interesting to investigate the stability properties of the solution (pq) . To this purpose we 
consider deviations from this stationary value: 



PQ 



PM = /PM ] +6 



PQ 



(87) 



Since {pm/pq) b = 0, the resulting equation is 

S = 3H 
By visualizing this dynamics in a 5 — 6 diagram, the stationary solution (|86j) is an attractor solution for 



Pm 
PQ 



Pq_Pm 
Pq Pm 



P J1 

PQ 



P 



M 



^- — <0 

PM 



r N 

o<- <r = — 

r N 



(89) 



The rate N/N represents a limit for the rate r/r at which energy is transfered from the Q component to the matter. 
We conclude that a certain decay of vacuum like Q matter into CDM offers a potentially promising approach to 
the coincidence problem. The phenomenological concept of a decaying vacuum is widely used in the physics of the 
early universe in order to describe the decreasing dynamical role of an effective cosmological "constant" connected 
with a transition from an initial inflationary period to a subsequent FLRW behavior [ pO| , pl[ . It is note- worthy that 
a similar mechanism seems to be relevant also for the transition from matter dominance to "vacuum" dominance in 
the late universe. At first sight it might seem counter-intuitive that a decay of the vacuum at the same time leads to 
an apparent dominance of the latter. However, for the previously discussed model B oc —H^ 1 (Sec. IV. C.) we shall 
confirm explicitly that such kind of transition is indeed consistent with a positive rate r/r > 0. 

At this point we come back to our previous statements (following Eqs. (g) and (jig)) on the nature of "particle 
production" in the presented formalism. According to Eq. ( |75| ) the ratio \B\/H, equivalent to T/3H determines the 
contribution of the "vacuum" to the total energy density. Cosmic anti-friction may be viewed as a vacuum effect 
which is connected with a nonvanishing "particle production" rate T. For an equation of state pm = Hmw it is 
obvious, that r/r is the production rate of matter particles out of the decaying "vacuum" component Q. It represents 
a real "physical" particle production. According to Eq. (p9T) the rate r/r is smaller than T = N/N. Only for 
the "stationary" solution (pq ) the "physical" particle production rate coincides with T. For r — const the "particle 
production" is entirely connected with the Q component. 

In our basic setting the cosmic substratum is entirely made of non-relativistic particles which are governed by 
an equilibrium distribution function. The splitting ( f73[ ) has revealed that this "generalized" equilibrium description 
may be regarded as a two-fluid model of "conventional" and "exotic" matter. The point is that also the "exotic" 
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matter is described in terms of "conventional" matter particles, only that the property of being "exotic" requires a 
"production" process. If generated at a certain rate, "conventional" particles in the expanding universe effectively 
exhibit vacuum-like properties, which manifest themselves macroscopically through a negative pressure. 

Our fluid approach formally implies that a particle number is attributed to the "vacuum" component as well. At 
first sight the concept of a particle number of the vacuum might appear obscure. In a fluid picture, however, it appears 
quite naturally [cf. Ref. |5jJ ], as long as the interpretation of the first moment (14) of the distribution function as 
"particle number flow" is maintained. As in any two-fluid picture, the production rate T naturally splits into 



nT = njv/r 



M 



l Q l Q 



(90) 



implying a split of the total particle number density n of the one-component model into n = nu + "q, by which 
the notions of a particle number density uq of Q "particles" and a corresponding change rate Tq are introduced. 
Tm = r/r is the rate by which the CDM particle number changes, Um is the CDM particle number density. A 
non- vanishing T does not necessarily imply a non- vanishing T m but may be a feature of the fluid picture of "vacuum 
matter" , characterized by Tq ^ 0. 

B. The case B oc -H' 1 

In this subsection we focus on the previously discussed case B ex —H^ 1 [cf. Eqs. (|66|)-(f7l|)] for which the energy 
density is given by 



Po 



p(l + zf + l 



(91) 



Performing the splitting (|73|), it is convenient to replace r by 



r{z) = 



[X + l 



/(*) 



'"o = 



H + l 



The energy density (£H|) then decomposes into 



Pm 



H+l 



Po (l + zTf(z) 



PQ 



Po 



/o = l 



l-(/(z)-l)/i(l + *)' 



(92) 



(93) 



For the special case / = 1 this splitting characterizes a non-interacting mixture of non-relativistic matter and "vac- 
uum". This is just the ACDM model as already mentioned below Eq. (f39|). In the general case / = f (z) both 
components interact. More specifically, there will be a decay of the Q component into matter. The constant p is the 
ratio 



t< 



p M (z = 0) 
PQ (z = 0) 



(94) 



Here it is expedient to emphasize that the basic one-component dynamics ( p7[ ) and (UW is compatible with any 
continuous positive-definite but otherwise arbitrary function /. The splitting into two components does not affect 
the quantity II/p at all. However, different splittings may produce different perturbation spectra. In particular, 
isocurvature perturbations may occur which should be sensitive to the type of splitting. Corresponding effects are 
expected to leave an imprint on the CMB anisotropics and are potential tools to discriminate between different choices. 
Furthermore, the underlying two-component dynamics is relevant for the coincidence problem. The corresponding 
stationary solution pm/pq = const is easily found. The condition (JS5|), together with Eqs. ( P2[ ) and (|7l|), provides 



/.= 



1 



/'- 



1 



H + (l + z)" 



(95) 



where the subscript s again indicates stationarity according to Eq. (pq). The obtained solution (p3) for / is expected 
to describe the splitting (B3) for large cosmological times. For the corresponding asymptotic behavior of Pm and pq 
we obtain 



(Pm) s 



PPQ 

iP+lf 



1 + n (1 + z)" 



(96) 
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and 



A) 



v! 



G°Q) 8 = r^5 l + M(l + *r > (97) 

(M + 1 ) l j 

respectively. These asymptotic expressions satisfy 

^-) = ix = -^- = const . (98) 

PQ / s 1 ~ r o 

Both components red-shift at the same rate. There is a permanent energy transfer from the Q component to the 
matter component. Without transfer pm would red-shift as a~ 3 while pq would remain constant. The transfer makes 
Pq red-shift also and, on the other hand, pu to red-shift at a lower rate than without transfer. In general, the red- 
shifts differ. For a specific amount of transfer, however, given by the expression ( Bq ) for /, the rates just coincide. We 
conclude that a fixed ratio pm/ PQ is compatible with cosmic anti- friction. If our universe follows such an asymptotic 
solution presently, this should manifest itself by a violation of the CDM particle number. If we are still far from 
the asymptotic regime, the CDM particle production might be negligible. These considerations of the coincidence 
problem give rise to a picture according to which an initial cosmological term did not completely vanish during the 
early stages of the cosmological evolution but still exists and even continues to decay. At large cosmological times not 
only the cosmological term itself becomes dynamically relevant again but also its decay properties may be essential 
for the asymptotic state of the universe (ignoring here a recently discussed scenario with a less bleak eschatological 
picture |p2| ). 

VI. SUMMARY 

We have introduced the concept of cosmic anti-friction to discuss the possible origin of an accelerated expansion 
of the present universe. Cosmic anti-friction relies on a one-component picture of the cosmic substratum which 
is regarded as a self-interacting gas of non-relativistic particles in (generalized) equilibrium. Together with simple 
assumptions about the interaction rate it allowed us to establish exactly solvable models of the cosmological dynamics. 
Cosmic anti-friction leads to a negative bulk pressure which may well account for the magnitude-redshift data of type 
la supernovae. The ACDM scenario is recovered as a special case of cosmic anti-friction. For the models of subsections 
IV. B. and IV. C, which are consistent both with the SN la observations and with CMB anisotropy data, we find the 
beginning of the phase of accelerated expansion at redshifts z acc rs 1.8 and z acc « 0.7, respectively. The one-component 
dynamics may be split into a two-fluid mixture in which a Q matter component decays into CDM. For a suitable decay 
rate there exists an attractor solution characterized by a fixed ratio of the energy densities of both components, which 
indicates a possible solution of the coincidence problem. As for scenarios with a cosmological constant or quintessence, 
the microphysical evidence for our models remains open. This shortcoming seems presently unavoidable and reflects 
our basic ignorance concerning the substance our Universe is made of. 
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